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Abstract. 

In present work the effective singlet-triplet model for Cu02-layer in the framework 
of multiband p-d model of strongly correlated electrons is obtained. The resulting 
Hamiltonian has a form of generalized singlet-triplet t-t'-J model for p-type 
superconductors and form of usual t-t'-J model for n-type superconductors. In the 
mean field approximation in X-operator representation we derived equations for Gorkov 
type Green functions. The symmetry classification of the superconducting order 
parameter in the case of tetragonal lattice resulted in d x 2_ y 2- and d xy -types of singlet 
pairing for both p- and n-type superconductors while s-type singlet pairing don't 
take place. Also normal paramagnetic phase of effective singlet-triplet model was 
investigated and the Fermi-type quasiparticle dispersion over Brillouin zone, density 
of states and evolution of Fermi level with doping were obtained. 
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1. Introduction 



Almost twenty years passed from the discovery of High-T c Superconductivity (HTSC) 
but for now there is no widely accepted theory of pairing mechanism in cuprates of p- 
and n-type. The necessary ingredient for discussing possible mechanisms of HTSC is 
the band structure of the fermion-like quasiparticles. However, it is a difficult subject 
for ab initio calculations due to the strong electron correlations. For this reason we will 
use a model approach and in order to get an agreement with experimental data we will 
start with a realistic multiband p-d model of transition metal oxides. 

The Hubbard model is often used to study electronic structure of strongly correlated 
electron systems (SCES). To take into account the chemistry of metal oxides the 
Hubbard model is generalized to the p-d model, a simplest version of such a model has 
been proposed by Emery |l| and Varma et al. ||. In this 3-band p-d model only d x 2_ y 2 
Cu and p a O orbitals are considered. One of the important features omitted in this model 
is asymmetry of n-type (electrons doped) and p-type (holes doped) cuprates. Point is 
that the spin-exciton concerned with singlet-triplet excitation of two-hole term occurs 
only in p-type systems, but not in n-type ||. Other feature is a non-zero occupancy 
of d z 2 Cu orbitals j|. There also a dependence between T c and occupancy of d z i was 
found. Hence more realistic model of CliCVlayer has to include d x 2_ y 2- and d z i- orbitals 
on copper and p x -, p y -, p z - orbitals on oxygen as well. Such model is the multiband p-d 
model that has been proposed by Gaididei and Loktev 
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Here r and i are Cu and O sites, A = d x 2_ y 2,d z 2 and a = p x ,P y ,Pz are orbital indexes 
on given copper and oxygen site respectively, e d and e p are energies of d x 2_ y 2- and d z 2- 
holes on copper and p x -, p y -, p z - holes on oxygen. U d , U p are intra-atomic Coulomb 
interactions, t pd is hopping integral for nearest neighbors copper-oxygen, t pp is hopping 
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integral for oxygen-oxygen, V dd , V pp , V pd are inter-atomic Coulomb interactions and 
J dd , J pp are exchange integrals. Abbreviation "H.c." means Hermitian Conjugation and 
" < r,i >" denotes that sum runs only over indices r ^ i. 

Simplest calculations of this model have been made by exact diagonalization of 
CuOi || and CuOq clusters. It was shown that when c/^-orbitals are neglected then 
the triplet with energy e 2 T lies above singlet with energy e 2 s on about 2 eV. This lets us 
to neglect (i 2 2-orbitals and we immediately come back to the 3-band model. However, at 
approach of <i 2 2-orbital energy to rf x 2_ y 2-orbital energy the difference e 2 t — £2S decreases 
and, at the certain parameters, crossover of singlet and triplet take place. Similar 
results where obtained in Q by self-consistent field method and in || by perturbation 
theory. All these facts give us cause for deeper investigation of processes concerned with 
two-particle triplet. 



2. Formulation of Effective Model 



In this paper we will use Hubbard operators (or so-called X-operators) Xj q = \p) (q\ 
on site f instead of annihilation and creation operators because they are a good tool 
in the case of strong electron correlation. Also defining that m <-» (p, q) numerate a 
quasiparticle described by the Hubbard operator X pq , and ^\ a {jn) is a parameter of 
X-operator representation for the single-electron annihilation operator with orbital A 
and spin a we can have the next correspondence between annihilation (and creation 
operators) and Hubbard operators: 

a /ACT = $>Aa(m)Xf . (6) 

m 

Hermitian conjugation of Hubbard operator indicated by cross: X™ + . 

To calculate quasiparticle band structure in SCES the Generalized Tight-Binding 
(GTB) method has been proposed 0. This method combines the exact diagonalization 
of a cell part of the Hamiltonian and the perturbation treatment of the intercell part 
in the X-operator representation. In paper [|l(J the consequent development of GTB 
method for La 2 CuO^ with a CuOq cluster as elementary cell was given. The problem 
of non-orthogonality of nearest clusters' molecular orbitals was solved in direct way 
by constructing explicitly Vanier functions on d x 2_ y 2,d 3z 2_ r 2,p x ,p y ,p z - five orbitals' 
initial basis of atomic states. In a new symmetric basis one-cell part of Hamiltonian 
are factorized allowing to symmety classification of all possible effective one-particle 
excitations in Cu0 2 plane as transitions from n-th hole term to (n+l)-hole term. The 
X-operators are constructed in the Hilbert space that consists of a vacuum = state 
|0), single-hole \a) = {| T)> I !)} molecular orbital of b\ g symmetry, two- hole singlet 
\S) = | |) I) °f 1 Ai g symmetry and two-hole triplet \TM) (where M = +1,0,-1) of 
3 B\ g symmetry states. In the X-operator basis the multiband p-d model Hamiltonian 
is given by: 



H = T,Uj2 X 7 + ^sXf + e 2T £ Xj MTM ) 
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+ E V, 

Here the energies £1, e 2 g and e 2 r are counted off from chemical potential /x, and t/ g 
indexes 0, a and 6 relevant to quasiparticle in lower singlet (0), in higher singlet (6) 
and in higher triplet (a) Hubbard's bands. The scheme of the levels and available 
quasiparticle excitations between them are presented in figure |l| 




Figure 1. Quasiparticle excitations between state in full singlet-triplet model 
(excitations for only one spin projections are shown). 

The t-t'-J model is derived by exclusion the intersubband hopping between low 
(LHB) and upper (UHB) Hubbard subbands for the Hubbard model [1 (| and for 3-band 



p-d model ||12|| . We write the Hamiltonian in the form 

H = H + H U (8) 

where the excitations via the charge transfer gap are included in Hi. Then we define 
an operator 

H(e)=H + eHi (9) 
and make the unitary transformation 

H (e) = exp (~ieS) H (e) exp (ieS) . (10) 
Vanishing linear in e component of H (e) gives the equation for matrix S: 



Hi + i H , S 



0. 



(11) 



The effective Hamiltonian is obtained in second order in e and at e = 1 is given by: 



H = H + -i [Hi,S 



(12) 
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It is convenient to express the matrix S in terms of X-operators [ D| . Thus, for the 
multiband p-d model ([!]) with electron doping (n-type systems) we obtain the usual t-J 
model: 



-Uiflj i , 



(13) 



H t _j = ^eiX° a + E t0 ij x i° x j a + 2 ^ J V (s^i ~~ 

i,a <i,j>,a <M> 

where spin operators Sf and number of particles operators rif can be expressed in terms 
of Hubbard operators as follows: 

1 
4 



s,s, -r^n 3 = W [xfxf - xrxf) . 



(14) 



For p-type systems effective Hamiltonian has the form of a singlet-triplet t-t'-J 



model flTB] : 



H = H + H t + Hj 



(15) 



where H (unperturbated part of the Hamiltonian), H t (kinetic part of H) and Hj 
(exchange part of H) are given by the following expressions: 



E, 

and 5Jij is the correction to it (dependent on the triplet's contribution 



(16) 



= E U E + ezsX? s + e 2T E Xj MTM ) , 

i V <t M J 

H t = E t^Xf'Xf + E tf g (aV2X™ * - X^ °) (aV2XJ T0 - XJ T2 ° 

<i,j>,a <i,j>,a 

+ E tp^b [Xf ° (crV2XJ T0 - XJ T2ct ) + H.c] , (17) 

<i,j>,(T 

Hj = \Y. + (SiSj - \rnnj) -\ E SJ^Xf. (18) 

<i,j> V 7 <*J>,o- 

The Jjj is the exchange integral: 



Jii = 4-^-, (19) 



5 j.. = 2t; 2 ^-^. (20) 



,2 

For the nearest neighbors (i = 0, j — 1) the estimation gives: 

^ ~ 1(T 2 . (21) 

Here _E ci is the charge-transfer energy gap (similar to U in the Hubbard model, 
E ct « 2eV for cuprates). 

Previously the motion of triplet holes and the simplest version of singlet-triplet 
model have been studied in |b|, [TJ|]. But these authors used the 3-band p-d model 
|I|, |2j where the difference between energy of singlet and energy of triplet Et — £s is 
about 2 eV. This fact leads to negligible contribution of singlet-triplet excitations in 
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low-energy physics. But in multiband p-d model St — Ss depends strongly on various 
model parameters, particularly on distance of apical oxygen (p 2 -orbitals) from planar 
oxygen (p x - and p y -orbitals), energy of apical oxygen, difference between energy of d z i- 
orbitals and d^.^-orbitals {e& 2 — 2 _ 2 ). For the realistic values of model parameters 
st — Ss appears to be less or equal to 0.5 eV (see || and JT0|). It is this small value 
of singlet-triplet splitting that let us believe that singlet-triplet excitations has a non- 
negligible contribution to band structure and also can give new superconducting pairing 
channel. 

As easily can be noted the resulting Hamiltonian (fl5|) is the generalization of t-J 
model to account of two-particle triplet state. But the inclusion of this triplet leads 
to such significant changes in Hamiltonian as renormalization of exchange integral (^0|) 
and appearing of Xf cr XJ fJ term in Hj. 

More significant feature of effective singlet-triplet model (|15"D is the asymmetry 
for n- and p-type systems. This fact is known experimentally. In particular, the 
holes suppress antiferromagnetism more strongly than electrons. It was observed in 
La2- x Sr x Cu0 4: in compare with Nd2- x Ce x Cu0 4: |16j. For n-type systems the usual t-J 
model takes place while for p-type superconductors with complicated structure on the 
top of the valence band the singlet-triplet transitions plays an important role. In first 
case we have spin-fluctuation pairing mechanism (see review [17]]). In the second case 
in addition to spin-fluctuations described by Hj we also have pairing mechanism due 
to singlet-triplet excitations. Really, let's look at structure of H t . There we can see 
terms like Xf a X° TQ which can be identically written according to multiplication rule 
of Hubbard operators: 

X^XJ T0 = X^XJ S X? T0 . (22) 

Here we can see three processes: creation (Xf a ) and annihilation (XJ S ) of the hole 
at different sites i and j, and spin-exciton Xj T0 . This spin-exciton can play a role of 
intermediate boson in superconducting pairing. 



3. Green Functions for Effective Singlet-Triplet Model 

Processes described by Hamiltonian ([15]) are shown in the figure |2|. We have three 
Fermi-type excitations and hence three base root vectors cti, ct2, a^. Corresponding 
basis of the Hubbard operators is: 

{X! S ,X* T0 ,Xf =► {XIX?,X?} . (23) 

The condition of basis completeness has a form: 

*7 + xf + J2 X J M ™ = i, (24) 

cr M 

where M = {2a, 0, 2a}. 

In the introduced notations H t has the following form: 

Ht= E E7^K«)^- (25) 

<i,j>,a m,n 
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Figure 2. Quasiparticle excitations between state in effective singlet-triplet model 
(excitations for only one spin projections are shown). 



where m and n are root vectors. Also, for further convenience, direct dependence on 
spin cr is introduced. Matrix 7y j0 .(m, n) can easily be obtained from Hamiltonian (|T5|) : 



+bb 1 -tab 

1 +aa 



ab 



\ -2a lb t ab 



i.i 



L -t a 
-aV2t% 



-i°ibtt \ 



4-CLCL 



(26) 



Equations of motion for Hamiltonian ( JT5I ) can be expressed as: 
d 



i— X p 

dt fa 



X p faJ H + H t + Hj 



(27) 



and ftp = 



where L% = L%(t) + L%(J), L%(t) [X%,H t ], L%{J) = [X%,H 
[X P f a ,Ho\ is vector of one-electronic energies: 

' S 2S ~ £i ^ 

V p = e 2T -s 1 . (28) 

\ S 2 T - Si J 

Taking into account equation (|||) and identity [A, BC] = {A, B} C — B {A, C} we 
can right away get an expression for L p a {t): 



L %(t) =EE^K") (Ef(p,m)X^ -X™, + Df(p, 



n 



(29) 



where E a /{p,m) = {x p fa , X™, + } is the neutral boson and DJ a '(p,n) = [X p f(T , XJ a ,\ is 
2e charged boson: 



Ef (p,m)=5 a 



(xl 



SS + X aa 

xJ os 
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A / A f 







(30) 



/ 



D°f(p,n) = 0. (31) 

Vanishing of Dj a '(p,n) is due to neglect of two-hole excitations |0) — > \S) in the t- 
t'-J model. Nevertheless such excitations had their effect resulting in the superexchange 
interaction J. 

It's easy to rewrite expression for L P j- a (t) in reciprocal space (or so-called k-space): 



LiM = E E n)E^ q (p, m)X n q<T , 



(32) 



q,& 



Introducing for convenience operator: 



y t y G(J i X 7 y a(T 



to 
-q 



(33) 

(34) 
(35) 
(36) 

For decoupling of equations on Green functions we will use the method of irreducible 
operators |18|]. This method based on linearization of equations of motion within 
Generalized Hartree-Fock Approximation (GHFA). Let us first introduce the irreducible 
operator: 

_ - ' - ' ' (37) 



we can now write down L p a {J)\ 

rl t t\ \ ( yd~S y I v T ( V u ^ V^" _i_ y aa Y a S W 

^ka\ J ) ~ 2 \^k-q X Qar + 1 q<r J ^k-q ~ 2 <* \ k ~1 1 + 1 *k-q) ) 

T 2 ( i\ ^ \ "> ( ya TOy i y ya TO ^ T / x^cr TO yoo i ycra ya 

^ka\ J ) — 2 £^ y k-q i a i a k-q — 2 9 V 9 +A d A fc- 

L3 / T\ ^ \ s x^tj T2(Tax i V V " T2cr "^t / ya T2a yaa i x^cro' x^cr T2cr 
fco-Wj - p I A fc-g + Iq&^k-q -y J -q\^k-q A g + A g A fc- g 



where 



LI* = L p k . - E Cto(p, - E ^ 

C fc(7 (p, /i) 
A fcff (p, h) 



ka ' 



\L\(ji ^-ka] 



(3? 



(39) 



Here Afc CT connected to the superconducting order parameter. Its symmetrical 
properties will be analyzed in next chapter. 

Now we can write down expressions for Ck a {p, h) and A^Qo, h) in compact form: 

CUP, h) = ,^ n \, ^ E h)) 



A fc(T (p, h) 



7ff -fc,a(p, (E™(p,p)) + <f M (K_ 9CT (p, /i)) 
1 



(40) 



<££o(M)> 
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h-qA m > n ) Ri gJ (P> m 'i h > n ) -7- 3 #(™,n)#+ (2) (p,m;/i,n) 



(41) 



Tensors R$(p, m; h, n) and R q(T ^(p, m; h, n) are presented in [Appendix A 



Averaging-out in C ka (p, h) were made in Hubbard-I approximation and all averages 
in C ka (p, h) are as follows: 



SS I v 
k + A / 



V 



o 
o 





X fc T0T0 + Xf 




yT2crT2, 








(Y g a( P ,h)) 



( (Y q a) 

(Y q(T ) 
V (Y g ,) J 



(42) 



(43) 



where Y qcr is determined by ([331). Also the anomalous averages were introduced: 

B ka (p, h) = (x p _ k ,Xl) . (44) 

Direct calculations can reveal one useful property of this averages: 

B-k&(h,p) = -B kf7 (p,h). (45) 

At calculation of ({L 3 ka (J), X^S) the terms of type (x%J q 2 ° P^ q {hf} appears 
where P kq {h) = \fN[X q a , Xz kd ] is a Fermi-like operator. These terms are responsible 
for 3-fermion excitations and in the Hubbard-I approximation they are equal to zero. 

Tensors R q l }(p, m; h, n) and R q } 2 \p,m] h,n) are presented in Appendix I. Shortly, 
each of their elements is equal to anomalous average B qa or zero. Matrixes A^}_ k {p, h) 
and A q 2 + k (p, h) consist of exchange integrals and has the following explicit form: 



A q l l k (p,h) 



A%{p, h) 



( J„ 



\ 



\ 



(46) 



(47) 



lq-k Jq-k 2SJq-k 

Jq—k Jq—k 25Jq— k 

y 25J q - k 25J q ~ k 25J q - k J 

( Jq+k Jq+k ^ 
Jq+k Jq+k 

J 

We can now perform elementary check of the obtained A ka (p, h) and C k(T (p, h) by 
using fact that when moving energy of the triplet TM to the infinity we have to obtain 
the usual t-J model |TJ, 2(J. Because in paper |TJ| Hubbard operators were also used 
we will compare our gap and renormalization of spectrum with its results: 

J 



A k~a J 



Tin 



-rin 



fag ~~ o (J k +i + 



Tjt-J 



0, and t 1 }™ 



Neglecting all excitations to the triplet zone (ct2 — > 0, 0:3 — > 0, 5J g 
n or m is equal to a) coefficients and matrixes in singlet-triplet model becomes: 

7 /9j(T (m,n) = 5 ml S nl tf g , 



(48) 
(49) 
>0 if 

(50) 
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Having made all these limiting transitions in ([40D 
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(51) 
(52) 

(53) 



and (Ell) we derive: 



t b k \l-n 



(54) 
(55) 



This is exactly what we have in the t-J model. So the effective singlet-triplet model in 
the low energy limit (triplet energy tend to infinity) reduces to the t-J model. 

Let's return to decoupling of Green functions. In GHFA we have to put irreducible 
operator I7 ka to zero. In this case the equations of motion take the following form: 



i-Xl = £ [Q p S p>h + C ka (p, h)) Xl + £ A ka (p, h)X h _ k , 
at h h 

• dp + _ 
l dt X - k * - 



E [n p s p>h + c± ka ( P , h)} x\+ - y: ai w (p, h)x-, 



(56) 



h 

ka 



Now we can easily write down a system of Gorkov type equations for normal and 
abnormal, and this system is closed: 

E P " V P ) 5 p , h - CUP, h)} G ka (h, l)-J2 A fcCT (p, h)F£(h, i) = ({xl, xl + }) 
h r * . (57) 

E + <W + c^(p, /») + E JOGWfo = o 

Actually this is an equation for a matrix normal Green function G ka - 



vh 

A fc(7 



and matrix abnormal Green function F ka 



X 



h + 

-ka 



X, 



I + 

ka 



X, 



I + 



Introducing matrixes: 

= (Qp5 P ,h + C ka (p, h)) , 
A ka = A ka (p, h), 

and solving equations we obtain expressions for G ka and F ka . 



r ka 

G ka 



ka 



-1 



A± ka G ka 



Ei - §t ka + A ka (ei + 1 At J 1 



(58) 
(59) 
(60) 



(61) 



where I is identity matrix and N is number of vectors in k-space. 

It can be seen now by analogy with the BCS theory of low-temperature 
superconductivity that A ka is the superconducting order parameter. 
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The system fl6"T|) is the set of matrix Green function and can be used to obtain 
energy spectrum and averages for any problem with defined basis of root vectors. In 
following chapters we will use system (|Bl"D to perform symmetry classification of order 
parameter and to investigate normal paramagnetic phase. 

4. Symmetry Classification of Superconducting Order Parameter 

An early suggestion that AF spin fluctuation could give rise to singlet d x 2_ y 2- wave 
pairing in p-type cuprate superconductors was made by Bickers, Scalapino and Scalettar 



21]. This suggestion has been supported by the FLEX approximation to the Hubbard 



model |22] that is, however not valid in the case of U 3> t. In this limit of SCES the 



proper model is the t-J model. Exact diagonalization and quantum Monte-Carlo method 
results for small clusters have been discussed by Dagotto [p3]| . For the infinite lattice 
the most adequate perturbation approach to the t-J model has been formulated in the 
X-operator representation because of the exact treatment of local constraint due to X- 



operators algebra. The mean-field solution f24| , |25| , |26| of the t-J model and analysis of 
the self-energy correlations beyond the mean-field approximation by diagram technique 
p7| and by high-order decoupling scheme []20| has confirmed the d x i_ y i- pairing in the 
hole-doped system with typical T c [x) dependence. Latest experimental results appear to 
prove d x 2_ y 2- pairing not only in p-type systems but also in n-type systems (see review 

HI)- 



Let's proceed to a symmetry classification of the order parameter A ka (p, h) in the 
effective singlet-triplet model considering case of square lattice. First, we have to break 
hopping and exchange integrals in two terms: 

7a(Pj h ) = t(P, h ) u k + ^(Pj fy&k, (62) 

J k = Ju k + J'u k , (63) 

5J k = 5Ju k + 5J'£j k , (64) 

where 

oj k = ex P — 2 (cos k x + cos k y ) , (65) 

8 

u k = ex P (ik8') = 4 cos k x cos k y , (66) 

<5' 

and non-primed values are concerned with nearest neighbor (figure 0(a), first 
coordination sphere) and primed values are concerned with next-nearest neighbor (figure 
|3](b), second coordination sphere). 

Second, to make a classification we will distinguish following symmetry types: 

4-1. s-type 

In superconducting phase we have a constraint condition: 
1 

j~ r ( a k*a-ka) = 0. (67) 

^ k 
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(a) 



(b) 



Figure 3. Location of nearest neighbors (a) and next-nearest neighbors (b) on square 
lattice. 



Right side of this identity in the effective singlet-triplet model is equal to zero due 
to choose of basis of root vectors a,, i.e. because of absence of the transitions from 
lower to higher Hubbard bands. Moreover, condition ( |67| ) is satisfied only in case of p- 
and d-pairing but not for symmetric s-pairing. Hence in the singlet-triplet model the 
symmetric s-type singlet pairing is absent. 



4-2. p-type 

Triplet pairing of p-type is impossible because for realization of this symmetry there 
must be a ferromagnetic interaction (see e.g. fT^fl ) but in case of singlet-triplet model 
we have only antiferromagnetic exchange. 



4-3. d-type 



Singlet d-type pairing is forbidden neither by the constraint condition (|67D nor by the 
type of interaction. 

First consider d x 2_ y 2- pairing type. In this case there is a restriction: 

cos q x B qa (p, h) = -J2 cos % B qa(p, h), (68) 
q q 
and, as a consequence, 

^sing^^O, (69) 



and 



^sin q y B q = 0. 



(70) 



This immediately leads to the significant simplification of superconducting gap 
A ka (p,h) = At x2 - y2 \p,h): 



2A ( ?* 2 -« 2 \p,h) 



N(E™ (h,h)) 



(cos k x — cos k y ). 



(71) 



[d 2 2) 

where A CT x ~ v (p, h) is the impulse-independent part of the gap: 



J {V 1 / + V 2 /) /2 JV 2 / 



(72) 
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K' h = E( cos ^ - cos q y )B ga (p, h). (73) 
i 

Id 2 2) 

Impulse-independent A CT x ~ y (p, h) includes exchange integrals only for nearest 
neighbors and therefore can exist only in boundaries of first coordination sphere. 

Now lets consider d xy - pairing type. The restriction made by this symmetry is as 
follows: 



Also 



cos q x cos q y B qcT (p, h) = 0. (74) 



sin q x cos q y B q = 0, (75) 



and 



Y cos q x sin q y B q = 0. (76) 
g 

The superconducting gap Afc CT (p, h) = A^ y \p, h) takes the form: 
\M( M 4A^)(p, h) . . 

AL„ y '(p,h) = = sin^sinL. 77 



where impulse- independent part of the gap A^ xy \p, h): 

( J'® 1 / J' ($i< 2 + $2,1) / 2 5 J'$3,l \ 

A (d»,)= J/ ($1,2 + $2,1) / 2 J'$2,2 ^J/^3,2 ; 

v 5 J'® 1 / 5J'<S> 2 / 5J'$ 3 / j 

$% h = J2 sin q x sin q y B qa (p, /i) . (79) 



(7? 



Coexistence of d x 2_ y 2- and rf^- pairing types is forbidden by the group theory: 
in the considered case of tetragonal lattice symmetry d x 2_ y 2- and d xy - types belongs 
to different irreducible representation (see review [EH]) so there must be concurrence 
between these types of pairing. 

Note that triplet channel results in 5 J in matrix ( [7^ ) and 5 J' in matrix (|78|). 
Obviously triplet channel gives the additional pairing. 

It is very interesting that d^-type includes exchange integrals only for next-nearest 
neighbors (see equation (|78|) ) and therefore can exist only in second coordination sphere. 
Hence in nearest neighbor approximation we can have only order parameter of d x 2_ y 2- 
symmetry. 

Actually equations (|7l|) and (|77|) are equations for order parameters \E f £' /l and 
$p,h reS p ec tively and should be solved self-consistently via Green function equations 
(|6T|). When this is done we can obtain energy spectrum and phase diagram T c versus 
concentration x. 
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4-4- Other symmetry types 

The other symmetry types are not realized due to absence of corresponding combinations 
of trigonometrical functions in expression for superconducting gap (fill). 

Combination of previously analyzed types such as s+d is impossible because in 
the considered case of tetragonal lattice symmetry s- and d- types belongs to different 
irreducible representation (see review |28[] and references from there). So the order 
parameter symmetry must be of d-type or, more precisely, to be a concurrence of d x 2_ y 2 
and d xy singlet pairing types. This conclusion is true not only for hole doped systems 
but also for electron doped systems because in the limit of infinite energy of singlet- 
triplet excitation (i.e. absence of excitations to the triplet states - this case corresponds 
to n-type cuprates) our equations take the form of t-J model ones and the gap symmetry 
remains. 



5. Normal Paramagnetic Phase 



In normal paramagnetic phase Afc CT = and Green function equations ( |61~D become 
essentially simpler: 



Gka = [EI — !R fcc 



Er)/VN 



(80) 

7 Z = o 

By solving these equations one can obtain energy spectrum and self-consistently 
find Fermi level. Also using following definition for the spectral density in terms of 
one-electron annihilation operators c\: 

A(k,a,E) = ~lm((ct\c ka )), (81) 

we can calculate density of states (DOS): 

N ( E ) = lj2 A (k,cr,E) = 1 £ (-^Im((c+, \c kaX ))) , (82) 

fccr k,a,X 

where a is spin and A is the orbital index. In terms of Hubbard operators the expression 
for DOS has the following explicit form: 



N(E) = ~ £ [{ll + it) Im ((Xl + \XI 

kcr 



+ {il + it+il) {^{{xi; 

= ^E[(7, 2 + 7, 2 ) ImGL(M) + 



+ 



lm((Xl~ 



+ (7 2 + 7a + 



-ImGL(2,2) + ImGL(3,3) 



(83) 



Here e and iVo(e) is the dispersion and DOS for non- interacting case, 7 m is the coefficient 
of transformation Ck a = J2lmXJI^, index "a" of Green function denotes that G% a is 
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advance Green function (while all over the paper we have used retarded Green functions). 
There is simple relation between this two types: Gk a (p, h) = —Gl a (h,p). 

The calculations described above where performed for High-T c superconductor 
La 2 -xSr x Cu0 4: with the following set of the model parameters (in eV): 

t pd = 1, e(d x 2_ y 2) = 0, e(d z2 ) = 2, e(p x ) = 1.5, e(p z ) = 0.45, 

t pp = 0.46, t' pp = 0A2, U d = 9,U p = 4, V pd = 1.5, J d = l. 



The same parameters where previously used in |T(J for undoped Lai- x Sr x CuO ±. 
The energy dispersion, obtained there, proved to be in good agreement with 
experimental ARPES data. 

Our results for slightly overdoped copper oxide (concentration of dopant x = 0.2) 
are presented in figure £|. 



Energy (electron representation) over Brilloin zone DOS(E) 
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Figure 4. Energy Dispersion (on the left) and Density of States (on the right) in 
normal paramagnetic phase of the effective singlet-triplet model for concentration of 
dopant x = 0.2. Energy of the singlet sub-band and its DOS are marked by dash- 
dotted line. Energy of the triplet sub-bands and corresponding DOS are marked by 
straight and dotted lines. Dashed line indicates location of the Fermi level. 
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Interesting case of coincidence of Fermi level and Van-Hove singularity is shown in 
figure |. 



Energy (electron representation) over Brilloin zone DOS(E) 
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Figure 5. Energy Dispersion (on the left) and Density of States (on the right) in 
normal paramagnetic phase of the effective singlet-triplet model for concentration of 
dopant x = 0.43. All definitions are the same as in figure |[ Distance between Van- 
Hove singularity and Fermi level is 0.074 eV 

It can be seen that singlet sub-band is very wide - almost 3 eV. It is a consequence of 
Hubbard I approximation in spectrum renormalization term Cka{p,h) (equation (f40|)). 
In more rigorous approximations with inclusion of spin correlators the energy bands 



should be more narrow (see |29j and |30|l ). 

At the concentration x = 0.43 Fermi level coincides with Van-Hove singularity in 
singlet sub-band. But if we will neglect hopping on second coordination sphere (nearest 
neighbor approximation) then we will get the flattening of dispersion in (tc, 0) — (0, it) 
direction and shift of Van-Hove singularities to higher energies. This will result in 
coincidence of Fermi level and singularity in singlet sub-band at concentration x = 0.33. 
It is typical value for the t-J model with nearest neighbor hopping t. It is also known 
that in t-t'-J model with the next-nearest neighbor hopping t' the same shift in energies 
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takes place so this phenomena is common both to t-t'-J model and the effective singlet- 
triplet model. Accepting the Van-Hove scenario of superconductivity, where the optimal 
doping correspond to coincident of Fermi level with Van-Hove singularity, we can clearly 
see that this should happened at x = 0.43. Meanwhile the experimentally obtained 
optimal doping value for Ld2- x Sr x CuOi is 0.18. So, either the Van-Hove scenario is not 
applicable or the model should be refined. This question could be answered only after 
the complete theory of superconductivity in the frame of the effective singlet-triplet 
model will be constructed. 

6. Conclusion 

In present work we obtained the effective Hamiltonian of the singlet-triplet model for 
copper oxides. Being the generalization of the t- J model to account of two-particle triplet 
state resulting Hamiltonian has several important features. At first, the X[ a Xj a term 
appear which can have a non-trivial contribution to superconducting pairing. Second, 
the account of a triplet leads to renormalization of exchange integral. And third, the 
singlet-triplet model is asymmetric for n- and p-type systems. For n-type systems the 
usual t-J model takes place while for p-type superconductors with complicated structure 
on the top of the valence band the singlet-triplet transitions plays an important role. 
The asymmetry of p- and n-type systems is known experimentally. 

The analysis of the possible processes in the effective singlet-triplet model shows 
that besides spin-fluctuation superconducting pairing mechanism typical for the t-J 
model we have pairing due to singlet-triplet transitions. These transitions induce spin- 
exciton which can play a role of intermediate boson in superconducting pairing. 

We have also performed a symmetry classification of superconducting order 
parameter. It was shown that in case of tetragonal lattice s-type singlet pairing doesn't 
take place while the d x 2_ y 2- and d xy -types of singlet pairing can exist. Moreover, there 
must be a concurrence between d x 2_ y 2- and d xy -types. At the same time, d^-type can 
exist only within the second coordination sphere and ^2 -type can exist only within 
the first. This fact lets us to take into account only d^^-type symmetry of order 
parameter in nearest-neighbors approximation. 

As concerns n-type cuprates, the gap symmetry was unclear for a long time. 
Recently, phase sensitive tunnel experiments by Tsuei and Kirtley p8| find an evidence 
for dominant d x 2_ y 2 symmetry in electron doped cuprates. This fact coincides with our 
results because in the limit of infinite energy of singlet-triplet excitation (i.e. absence of 
excitations to the triplet states - this case corresponds to n-type cuprates) our equations 
take the form of t-J model ones and the d^.^-gap symmetry remain. 

For normal paramagnetic phase we have obtained the energy dispersion over 
Brillouin zone and calculated density of states. Evolution of Fermi level with doping 
is also found. Both singlet and triplet excitations contributes to the density of states 
which leads to appearing of two Van-Hove singularities. At holes concentration x = 0.43 
the Fermi level crosses the first Van-Hove singularity corresponding to singlet sub-band. 
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And at holes concentration x = 0.64 the Fermi level crosses the second singularity 
corresponding to triplet sub-band. 
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Appendix A. Tensors R^J (p, m; h, n) and R+J- 2 \p, m; h, n) 

Tensors R^ 1 J(p,m; h,n) and R^ 2 \p,m; h,n) where introduced as follows: 

Vn([X\„ E^ q (p, m)]Xl) = 6 aa ,R$(p, m; h, n )+5^R ( S(p, m; fc, n).(A.l) 

The "+" symbol in R+^ (p, m; h, n) means that all averages B qa (p,h) in it has a 
interchanged indices p and h. 



Table Al. R\J (p,m; h,n) 







m=l 


m=2 


m=3 








h=l 


Bga(l,n) 










p= 


1 


h=2 





S 9CT (2,n) 











h=3 


-B qrT (3,n) 














h=l 


B qr7 (2,n) 










p= 


2 


h=2 





B qa (2,n) 











h=3 





-B qa (3,n) 











h=l 








—B qa (l, 


n) 


p= 


3 


h=2 








-B qa (2, 


n) 






h=3 














Table A2. R q J ' (p,m; h,n) 







m=l 


m=2 


m=3 






h=l 


-B qa (n, 1) 










p=l 


h=2 


B gff (n,2) 












h=3 








B qa (n, 


1) 




h=l 





-B qa (n,l) 







P=2 


h=2 





-B qa (n,2) 









h=3 








B qcr (n, 


2) 




h=l 


B qa (n,3) 










P=3 


h=2 





B qrJ (n,3) 









h=3 
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